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ABSTRACT 

We present a stellar population study of three Hll galaxies (Mrk 36, UM 
408, and UM 461) based on the analysis of new ground-based high resolution 
near-infrared J, H and Kp broad-band and Br7 narrow-band images obtained 
with Gemini/NIRI. We identify and determine relative ages and masses of the 
elementary star clusters and/or star cluster complexes of the starburst regions 
in each of these galaxies by comparing the colors with evolutionary synthesis 
models that include the contribution of stellar continuum, nebular continuum 
and emission lines. We found that the current star cluster formation efficiency 
in our sample of low luminosity Hll galaxies is ~10%. Therefore, most of the 
recent star formation is not in massive clusters. Our findings seem to indicate 
that the star formation mode in our sample of galaxies is clumpy, and that these 
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complexes are formed by a few massive star clusters with masses >1O^M0. The 
age distribution of these star cluster complexes shows that the current burst 
started recently and likely simultaneously over short time scales in their host 
galaxies, triggered by some internal mechanism. Finally, the fraction of the total 
cluster mass with respect to the low surface brightness (or host galaxy) mass, 
considering our complete range in ages, is less than 1%. 

Subject headings: galaxies: individual (Mrk 36, UM 408, UM 461) - galaxies: 
stellar content - galaxies: dwarf - galaxies: star clusters - infrared: galaxies 



1. Introduction 



Hll galaxies are dwarf galaxies undergoing an intense episode of star formation (SF) 
that dominates their total optical luminosity. Most Hll galaxies overlap in their observed 
properties, with Blue Compact Dwarf (BCD) galaxies; indeed, many are found in both 
samples. Typically, they have relatively small physical sizes (with a few kpc in size), and the 
starburst regions spatially cover the visible extent of the galaxy, making it difficult to access 
the presence of an underlying (or host galaxy) stellar population. Originally, their low heavy 
element abundance and the non-detection of an old population have given ri se to the question 
of wh ether they may be presently forming their first generation of stars ( ISargent fc Searle 
19701 ). Recent works, however, have shown that mos t Hll galaxies seem to present an 



underlying population from previous episod es of SF (e.g.. lTelles Ill995t iPapaderos et al.lll996 



Telles fc Terlevichlll997l : ICairos et al.ll2003l ). So, the integrated light observed in H ll galaxies 
is formed by the contribution of two main components, the young stellar population with 
ages of a few Myr to tens Myr and an underlying population of intermediate (with ages of 
hundreds Myr) to old stars (with ages > 1 Gyr). 



The structure of the Interstellar Medium (ISM) found in H ll galaxies ( iLagos et al.ll2007l . 
and references therein) have profound implications for the origin of the present starburst 
and the dominant large scale mode of SF. While many high luminosity Hll galaxies show 
evidence of morphol ogical disturbances that may be associated with interactions or ongoing 
minor mergers (e.g.. iTelles et al.lll997l : iBergvall fc Ostlin II2002I ). less luminous and compact 
H II galaxies defy o ur attempt to find morphological signatures of an external triggering agent 
( ITelles et al.l 119971 ). In fact, Hll galaxies seem to be typically isolated and not associated 
with giant galaxies (Telles & Terlevi ch 1995), and their clu stering properties seem to be 



similar to those of normal galaxies (ITelles &: Maddox 



many attem pts to search for H I companions (iTavloi 



companions (iNoeske et al.ll200ll : iPustilnik et al.ll200l[ ). In these studies, it was concluded 



20001) . There have been, however. 



19971 ) and intrinsically faint optical 
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that galaxy interactions with both massive and dwarfs are probably the main mechanism that 
triggers SF bursts in BCD progenitors. Although interactions are often invoked to explain 
burst of SF in Hll/BCD galaxies, it is possible that internal processes (e.g., gravitational 
cloud collapse and/or infall in conjunction with small perturbations) have been responsible 
for triggering the present episode of SF in at least a significant fraction of dwarf gal axies 



see 



Hunter fc Elmegreen 112004 : IPelupessv et al.ll2004j : lTellesll2009l : ISimpson et al.ll201ll and 



references therein). 

It has been apparent for over a decade now, with the advent of the HST, t hat the star- 
burst regions in these galaxies are composed of a myriad of star clusters (e.g., iBillett et al 
20021 ). with masses ^10^~^ MfD and sizes of a few pc, that are typically more massive than 



normal clusters in our Galaxy (IMeurer et al.lll995l : IConti &: Vaccalll994[ ). These Super Star 
Cluster (SSCs) or Young Massive Clusters (YMCs) were found originally in classical starburst 
galaxies and in ga laxies with evident sig ns of interaction or merger such as in the Antennae 
NGC 4038/ 4039 jwhitmore erall Il999l). dwarf irregular galaxies such as NGC 1569 and 



such as M82 flMelo et al. 



NGC 1705 (lO'Connell et al.lll994j) . and other sta r-forming dwarf galaxies or BCD galaxie s 



20051) and SBS 0335-052 (iThuan et al.lll997l : IPapaderos et al.lll998h . 



The formation of these massive clusters and the cluster mass function (CMF) is directly con- 
nected to the SF processes in the galaxies, in the sense that different physical factors, in these 
galaxies, put constrains on the mass of the star clusters. The CMF appears to be a single 
power law (with index ~-2), which implies the same SF mechanism for the massive clus- 
ters and their lower mass analogues. This type of young and massive c l usters were possibl y 



formed in high pressure conditions (e.g., lElmegreen fc Efremovl 119971 : iBillett et al.l |2002| ) 



hence the extremely high pressure regions give rise to more massive and compact clusters. 
We know that SF in clusters is a common phenomenon in starburst galaxies, and that the 
massive clusters play an important role on the evolution of the ISM of their host galaxies, 
producing large scale structures such as supershells (or bubbles) and creating galactic winds 
that cause, in some cases, the blowout of freshly produced metals from the galaxy into the 
intergalactic medium (IGM). The similar ity in mass a nd size between SSCs and Milky Way 
globular clusters (GCs; M~ 2 x lO^M©, lHarrislll99ll ) suggests the possibility of an evolu- 
tionary connection in the sense that SSCs are proto-globular cluster systems. The classic 
example of a YMC is R136 at the center of the 30 Doradus nebula in the Large Magellanic 
Cloud. This cluster has a diameter of ~1.6 pc, an age of 1-2 Myr (IMassey &: Hunten 119981 ) 
and an estimated mass of ~2xl0^ Mq (jWalborn et al.ll2002l . and references therein). 



Apart from evidence for massive star clusters and its effect on the surrounding medium, 
H II galaxies show an underlying low surface brightness (LSB) component likely to represent 
their host galaxies, whose stellar population is a product of previous episodes of SF. This 
underlying stellar host generally extends a few kpc from the central star-forming regions. 
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showing regular and elliptical isophotes with red colors indicative of an evolved stellar pop- 
ulation. The host galaxy in Hll galaxies has been st udied in the literature based on the 
analy sis of both optical (e.g.. iTelles fc TerleyichlllQQTl) and near infrared ( near-IR) images 



2000 



Cairos et all boosi iNoeske et al. l boosi boosi ) and spectroscopy (e.g., iRaimann et al.l 



Westera et al1bo04[ ). Ixhuanl dlQSsh observed, in the near-IR, a sample of BCD galax- 



ies and showed that the light in these galaxies is due to the presence of K- and M-giant stars. 
The study of the structural parameters of the host galaxy (exponential scale length (ao) and 
central surface brightness (/io)) has been commonly used in the literatur e to evaluate the 



evolutionary relations between the different types of dwarf galaxies (e.g., iPapaderos et al. 
19961 : iTelles et al.l 119971 : iGil de Paz &: Madorell2005l and references therein). The derivation 
of the ages and spatial distribution of these stellar populations is the first step towards 
establishing the evolutionary state and the SF history of these galaxies. 

We present near-IR broad-band J, H and Kp and Br7 narrow-band images of three low 
luminosity H ll galaxies: Mrk 36, UM 408, and UM 461. Our aim is to describe the properties 
of the star clusters or complexes which are distinguishable with our superb ground-based high 
spatial resolution images with NIRI on the Gemini North telescope as well to determine the 
structural properties of the underlying galaxy using surface photometry. These observations 
in combination with a proper assessment of recent stellar population synthesis models allow 
us to put some constraints on the recent and past history, and the dominant large scale mode 
of SF in these galaxies. 

The paper is arranged as follows: In § [2] we describe the sample, the observations, 
and data reduction. In § [3] we present the results and derive the properties of the star 
clusters/complexes and the LSB component. In §|l]we discuss our results and finally in §[5] 
we summarize our conclusions. 



2. Sample, observations and data reduction 
2.1. Our sample 

We targeted two galaxies Mrk 36 and UM 461, that are particularly rich in giant star- 
forming knots and one compact Hll galaxy, UM 408, with less evidence of multiple knots of 
SF. The light collecting power of the Gemini North telescope and the high spatial resolution 
permitted by NIRI are an excellent combination for resolving the star cluster or/and star 
cluster complex populations in these compact galaxies. In the following paragraphs we 
describe the main properties of our sample of galaxies. 

Mrk 36 (Haro 4, UGCA 225) is a compact Hll galaxy showing at least two large star- 
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forming knots (ILagos et al.l 120071 ) in the south-eastern part of the galaxy. iThuanI (119831 ) 
reports integrated ne a r-IR colors (de-reddened) of J — H = 0.38 and H — K = 0.42, while 
Hunter fc ElmegreenI (j2006[ ) reports integrated colors of J — if = 0.422 ± 0.100 and H — 
K = 0.543 =fc 0.310. The regio ns of largest H/3 line emission flLagos et al.l 120071 ) and HI 



maps (IBravo-Alfaro et al.l 120041 ) coincide, indicating that the current SF is restricted to the 
dense region of the parental cloud. Bravo-Alfaro et al. argue that a transient encounter 
between Mrk 36 and the neighboring spiral Haro 26 could explain both the SF in the former 
and the pronounced warp in the latter. Finally, r ecent radio observations with the VLA 
in 1.4, 4.9 and 8.4 GHz by iRosa- Gonzalez et al.l (120071 ) have shown that Mrk 36 has a 
ne arly flat radio spectruna domin ated by thermal emission, similar to the regions detected 
by iJohnson &: KobulnickyI ( 120031 ) in Henize 2-10. This may be an indication that the first 
contribution of the synchrotron emission to the low frequency emission, due to the first 
supernova (SN) explosions, has not yet appeared. These integrated properties show definite 
evidence of a very young starburst. 

UM 4O8 is a compact galaxy with a projected size of ^Ikpc. Although this galaxy was 



resolved as a single Hll region in previous studies (e.g., iGil de Paz et al. J2OO3I), two giant 



regions in the central part of the galaxy were detected by Lagos et al. ( 20091). with a ges o f 
~5 Myr and stellar masses of ^IO^Mq. Using GMOS-IFU observations iLagos et al.l (120091 ) 
showed that the metal content in the ISM of this galaxy is well mixed and homogeneously 
distributed throughout the galaxy, in the same way as in other dwarf galaxies. 

UM 461 (PGC 037102) is a well studied Hll galaxy. This object has a very compact an d 
bright off-center nucleus, some small regions spread along the galaxy (iNoeske et al.l 120031). 



and a n ext ernal envelope that i s strongly distorted towards the south-west. iDoublier et al. 
( I2OOII ) and iNoeske et al.l (120031 ) have obtained integrated near-IR colors that differ signif- 
icantly from each other, with J-H=0 .99 and H-K=-0.68 and J-H=0.47 mag and H-K=0.2 
mag, respectively. iTaylor et al.l (119951 ) proposed that UM 461 was formed together with UM 
462 by tidal interaction. But the HI maps of UM 461 and UM 462 do not show that these 
galaxies are tidally inter acting, therefore it i s unlikely that these objects induced the star 
formation on each other (jVan Zee et al.lll998l ). 



Finally, the presence of He ll A4686 emission bump of Wolf-Rayet (WR) stars has been 
reported in the integrated spectra of Mrk 36 and UM 461 (jContilll99ll ). 
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2.2. Observations, data reduction, and calibration 

Broad-band J(1.25/im), H(1.65/im), Kp(2.12/im), and narrow-band Br7(2.17/im) images 
were obtained using the NIRI instrument on tlie Gemini Nortli telescope on 2005 August 02 
(UM 408), November 24 (Mrk 36), and December 29 (UM 461). We used the f/6 camera 
which provides a field of view of -120"xl20"using the 1024x1024 pixels ALLADIN InSb 
detector, with a pixel scale of 0".116 on side. The observations were performed under pho- 
tometric conditions. Table [1] lists all observational parameters and values adopted in this 
work. In this Table, Column (1) gives the object name. Columns (2) and (3) give the a 
and 6 coordinates (J2000), respectively. Column (4) gives the observed heliocentric velocity 
(vel.) and the 3K CMB (Cosmic Microwave Background radiation) corrected distance from 
NED. Columns (5) and (6) give the oxygen abundance and extinction c(H/3) for these galax- 
ies obtained from the literature. Colurn n (7) gives the galactic extinction E(B-V) from the 



extinction maps of ISchlegel et al.l (119981 ). In Column (8) we present the date of observation. 
Column (9) shows the filter used in each observation and Columns (10) and (11) give the 
exposure time in seconds, considering the different coadd exposures and the mean air mass 
of each observation, respectively. Finally, Column (12) shows the instrumental constants Ca 
of the transformation equation for flux calibration as described below in this section. 

The data were reduced following the standard procedures for near-IR imaging using 
the Gemini/NIRI package version 1.8 inside IRAfQ. For each fllter, a normalized flat was 
constructed from the flat images observed with the calibration unit with the shutter closed 
(lamps off) and with the shutter open (lamps on). The bad pixel mask was constructed 
by identifying the bad pixels in the flat images with shutter off. The sky images were 
constructed from the raw science images by identifying all objects in each frame, masked 
out, and averaging the remaining good pixels (images were observed with a dither offsets 
of 10" and 20" for all galaxies). The raw science images were processed by subtraction 
the sky on a frame-by-frame basis and dividing by the normalized flat fleld images. The 
flnal flat-fielded, sky subtracted images were then registered to a common pixel position and 
median combined. 

In order to obtain a photometric calibration, we observed two standard stars: FS130 
(GSPC P264-F) and FS21 (GD140) for the observing run of Mrk 36, FS4 (SA93-317) and 
FS6 (Feige 22) for UM 408, and FS19 (G162-66, LTT 3870) and FS20 (G163-50, LTT 4099) 
for the observations of the galaxy UM 461. For each filter, the photometric calibration was 



^IRAF: the Image Reduction and Analysis Facility is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Universities for Research in Astronomy, In. (AURA) 
under cooperative agreement with the National Science Foundation (NSF). 
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obtained using the following equations: 

mxo = mx, + xK\ + Ca, (1) 

mx^ = -2.5 X log{c/t), (2) 

where m^o is the magnitude in the standard system, m^- is the instrumental magnitude 
with c the number of counts and t the exposure time, x the air mass, Ka the extinction 
coefficient and Ca the instrumental constant (A=J, H, Kp). The instrumental constant Ca 
was calculated as the average value from the standard stars. Given the limited number 
of standard stars observed, we have used the av erage value of extincti on to the "Mauna 



Kea" observatory (MKO) Kx=Q.l, 0.06 and 0.09 flKrisciunas et al.l 119871 ) to the filters J, H 
and K, respectively. The instrumental Ca constants used for each galaxy are listed in the 
last column of Table [H To compare the observed Kp measures with evolutionary synthesis 
models, we transforrned th e Kp into K magnitudes using the relation Kp= K + 0.22x(H-K) 



( iWainscoat fc Cowielll992l ). We selected star clusters/complexes by visual inspection of the 
images and using daofind in DAOPHOT. We fitted the detection threshold properly, with a 
value of 2.5cr above the local background, in order to detect sources in the starburst region 
that in the case of Hll/BCD galaxies spatially cover the visible extent of the galaxies. All 
the clusters or complexes in this study have been detected in the J, H and Kp band. The 
final catalog of objects is shown in § 13.11 



The Br7 images were calibrated using a procedure similar to that described in lLagos et al 



( 120071 ). We convolved the stellar Spectral Energy Distribution (SED) with the response of 
the narrow-band filter. To do this we used the SED, appropriately scaled, of stars with the 
same spectral type as our standard stars obtained from the hterature. Finally, we subtracted 
the continuum, estimating their contribution from the modeled SED for a given age of the 
regions obtained below in § 13.1.21 

Table [2] shows the observed integrated magnitudes for all the galaxies, uncorrected 
for Galactic extinction. The magnitudes were calculated by measuring the fiux inside fixed 
apertures of radii 14", 6", and 12" for the galaxies Mrk 36, UM 408, and UM 461, respectively. 
Foreground objects have been masked out. 
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2.3. Completeness limits 

We quantified our detection limits and test the reliability of the derived magnitudes 
using a series of completeness tests by adding artificial extended objects to our images. First, 
we created an empirical point spread function (PSF) for isolated point-like sources in the 
NIRI images. Given that the observed field of view (FOV) is relatively small (~120"xl20"), 
there is a possibility that these point-like sources are non-stellar. So these objects were 
selected by eye and we ruled out sources with irregular PSF shapes. Using this information 
we added 20-50 artificial extended objects distributed in a regular square array. Models of 



these extended sources were generated using mksynth in BAOLAB (ILarsen Ill999l ) and were 
added to the science images. Magnitudes were randomly assigned to each position from 16.00 
mag to 23.00 mag with an interval of 0.5 mag. Finally, we compared the number of added 
and recovered sources in each galaxy. The estimated completeness fraction for objects as a 
function of magnitude for the Kp band is shown in Figure [H The completeness goes down 
to 20.2 mag in Mrk 36, 21.5 mag in UM 408 and 22.03 mag in UM 461 with ~90 % of the 
objects recovered in the Kp band. 




Fig. 1. — Completeness profiles (Nrec/^add) in the Kp filter for star clusters in each of our 
studied galaxies, from the left to the right, Mrk 36, UM 408 and UM 461, respectively. The 
dashed lines show the recovered fraction of 90% and 50%. 



3. Results 

Figure [2] shows the three galaxies observed and studied in this work (Mrk 36, UM 408, 
and UM 461) in the Kp and Br7 bands, respectively. All galaxy images reveal the presence 
of bright regions and/or star cluster complexes, given the apparent sizes of these regions in 
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Figure [21 surrounded by LSB envelopes. Br7 emission is detected in all galaxies, but only the 
brightest regions show clear evidence of intense emission given by our detection limit. The 
unprecedented high spatial resolution images obtained for Mrk 36 allow us to identify, for 
the first time using ground-based telescopes, the elementary structures within the starburst. 
Mrk 36 shows a plethora of small clusters, distributed over the entire extension of the galaxy, 
but with a previously unresolved concentration in the central knot. The Br7 emission shows 
a peak in this region, indicating the young nature of these clusters. We also observed the 
presence of other groups of bright clusters in the Kp image located in the northern region of 
the galaxy with weaker Br7 emission. Here we named these two groups of regions as complex 
I and complex II, respectively. The galaxy UM 408 appears quite compact and regular, and 
the Kp-band image shows that the star-forming regions are distributed in a clearly defined 
ring-like geometry, composed of five almost regularly spaced regions where the current SF 
is occurring. The largest region is located in the eastern part of the galaxy. We observe the 
presence of another region outside the central part of the galaxy. The Br7 emission in this 
galaxy is very regular and weak. Fin ally, the Kp-band rn orphology of UM 461 is similar to 
that reported in previous works (e.g.. iNoeske et al.ll2003l ). The brightest region is off-center 
and located in the eastern part of the galaxy. This region is extremely compact and bright 
in Br7. We detected a number of regions spread throughout the body of the galaxy. 



3.1. Properties of the star clusters/complexes 
3.1.1. Near-IR colors 



For the color analysis we considered the objects obtained from our catalog in § [2] that 
are brighter than m^^-p ~21 mag. We detected 33 regions in Mrk 36, 6 in UM 408 and 13 in 
UM 461, respectively. All regions identified in the Kp-band images of the galaxies are shown 
in Figure El Given the seeing of ~0'.'4-0'.'5 we expect these regions to have diameters less than 
25 pc in Mrk 36, 88 pc in UM 408, and 37 pc in UM 461. So these diameters are upper limits, 
and the elementary cluster population in our samples of galaxies are barely resolved. Many 
of the detected regions are not individual star clusters but rather blends of several individual 
star clusters with s izes similar to that in other star cluster complexes in the literature (e.g., 
Melena et al.l l2009l ) . In the case of Mrk 36 the detected regions are compatib le with the 



sizes of some young star clusters detected in other Hll/BCD galaxies (e.g., M82; iMelo et al. 



20051 ). so these regions can be considered as individual star clusters, although likely some of 



— --^^ 1 "rt - 

( Zhang et al. 


2001; 


Larsen 


2004) 



derived are the luminosity-weighted mean value of the complex. We have to bear in mind 
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this caveat in our analysis but we are still able to derive the light weighted properties of 
these knots of SF and evaluate the dumpiness properties of t he large scale m ode of SF in 
these galaxies. In Figure [3] we also show the size (~200 pc; IWalbornl 1 199 if ) of the whole 
nebular region of 30 Doradus to compare with the sizes of the clusters in Mrk 36. 

We measured the flux of the individual star cluster /complexes in all filters (including 
Br7) using circular apertures with the program APER in IDL (an adapted version of the task 
DAOPHOT in IRAF). For each aperture, we considered only pixels with 3a flux above the 
background. For each cluster we calculated the colors J-H and H-K (after the transformation 
of Kp into K magnitudes). Figure H] shows the color-color diagram (J-H vs. H-K) for all star 
clusters/complexes. In this f igure we illustrate the evolutionary tracks of thes e colors using 



STAR BURST99 (black line; iLeitherer et al.lll999r ) and GALEV (orange line; iKotulla et al. 



20091 ) single stellar population (SSP) models, for a Kroupa IMF, and metallicity Z=0.004 
(more details in § I3.1.2p . The triangles represent the observed values for clusters in Mrk 36, 
the stars represent the values for the star cluster complexes in UM 408, and circles represent 
the values for the star cluster complexes in UM 461. Filled symbols indicate the detection of 
Br7 in these regions. Finally, we compare our observed apparent magnitudes with the ones 
found in the literatur e. The only o n e clus ter /complex found in the literature is the region 
#2 in UM 461 where iNoeske et al.l ( 120031 ) obtained mj=17.30 mag, that agree within the 
errors with our value of mj=17.34 mag. 

Table [3] lists the photometric values for each star cluster/complex identified in our 
sample of galaxies (see Figure [3]). In this Table, Column (2) shows the identification number. 
Column (3) shows the aperture considered to obtain the photometry in arcsec. Columns (4), 
(5), and (6) give the observed photometry in the J, H, and Kp bands for each cluster, 
respectively. The ^ symbol over the identification number indicates that Br7 emission was 
measured. Finally, Columns (7), (8), and (9) give the extinction E(B-V), age in units of 
Myr, and the stellar mass in units of IO^Mq of each star cluster or complex derived from 
the best fit comparison of the observed colors with chosen evolutionary synthesis models as 
described below in § 13.1.21 



3.1.2. Obtaining the physical properties of the star cluster/ complex population 

Models 

Absolute age dating of young star clusters in starburst galaxies has proven to be a very 
challenging task, due to the additional emission by dust plus gas, as well as the effects of 
extinction and metallicity. In relative terms the determination of ages and masses of the star 
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cluster population within these galaxies can show some insight on the overall dominant mode 
of recent SF. In order to show these effects we derived the extinction and ages using our 
color-color diagrams (see Figure H]) by comparing the near-IR colors for each star cluster or 
star cluster complex with two independent SSP models: STARBURTS99 models (model I; 
this model includes pure stellar and nebular continuum) and GALEV (model II; this model 
includes stellar continuum, nebular continuum, and emission lines from warm ionized gas) for 
an instantaneous burst of SF. So we derived the ages by calculating the best fit of the models 
(using the chi-square method) to the observed colors for the model, varying the extinction 
E(B-V) from to 1 in steps of 0.05 mag. We used the Galactic extinction curve given by 



Cardelli et al.l ( 119891 ) assuming Ry=3.1. The star cluster/complex masses were derived using 
the absolute magnitudes of each cluster as compared with the models. In the following 
paragraphs we describe the main properties of the models used in this work. 

Model I: We considered the STARBURST99 model, which includes pure stellar and neb- 
ular continuum, for metallicity Z=0.004 (^l/SZ©), and Geneva evolutionary stellar tracks, 
assuming a Kroupa IMF (ocM~") with a=1.3 for stellar masses between 0.1 to O.SM© and 
a=2.3 for masses between 0.5 and lOOMoJor^a total mass of IO^Mq. More details about the 



physics of the model in iLeitherer et al.l (119991 ). 



Model II: We compared our data with the GALEV evolutionary track for metallicity 
Z=0.004 and a Kroupa IMF (O.l-IOOM©) including stellar continuum, nebular continuum, 
and gas emission. The GALEV models used in this work were kindly provided to us by Ralf 
KotuUa. These models were run using the Geneva evolutionary tracks with a minimum age 
and time resolution of 0.1 Myr, unlike the models available on the GALEV website which use 
the Padova isochrones and have a minimum age and time step of 4 Myr. The models also use 
a fraction of visible mass, which is used for the standard models available on the Web page, 
and which are thus twice as bright as the Padova models for the same mass. This does not 
affect any of the colors, but we must keep this in mind when comparing masses. The flux of 
the hydrogen lines were computed using atomic physics and the production rate of ionizing 
photons, whereas non-hydrogen line strengths are computed using met allicity- depend e nt liii e 



ratios relative to H/3. More details about the input physics is given in lKotuUa et al.l (120091 ). 



Extinction, masses and cluster/ complex ages 

The two models (model I and II) displayed on the J-H versus H-K color-color diagram 
(Figure H]) show a different path for ages <6 Myr. The track for model I and the extinction 
vector, in the corner of this figure, are perpendicular to the model II. Therefore, estimates 
of extinction could be significantly different when we use the two cases. In Columns (7), 
(8), and (9) of Table [3] we show the E(B-V), age and mass of the star clusters/complexes. 
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using the two different models described above. In Figure [5] we show the distribution of E(B- 
V), ages, and masses for the star clusters/complexes detected in this work. In this Figure, 
the black distribution corresponds to the results obtained using STARBURST99 and the 
orange distribution to the those obtained using GALEV. From this figure it is clear that 
the extinction distribution is shifted to higher values when we use the model that include 
the contribution of nebular emission lines, meanwhile, the ages estimated using this model 
(model II) appear to be slightly younger and less massive than those obtained using model 
I. 

We can ask now, what is the most appropriate model to represent the ages and masses 
of our sample of star clusters/complexes? The SED of the observed star clusters/complexes 
in some Hll/BCD galaxies show a clear excess in the near-IR, that is reflected in red H-K 
colors with r espec t to th e models. T h is ob se rved red excess had been reported earlier by 
Vanzi et all J200ol. l2002h : Iflunt et all J2003h jjohnson et all tooi ). and more recently by 
Reines et all J2008h and lAdamo et al.N201ol. jgOllI ) in the star cluster population of other 
Hll/BCD galaxies. In fact lAdamq et al.l (12010 ) show ed t hat the I and H band s are signifi- 



Reines et al. 



mm and lAdamo et al.l (l2010f ) showed that 



cantly affected by a similar excess, 
this red excess clearly introduces a systematic offset in all of the derived parameters when all 
filters, from UV-optical to near-IR are included in the determination of the observed SEDs. 
This suggests that the models commonly used in the literature to calculate the properties 
of the c lusters, models that include pure stellar continuum, are inadequate for age-dat ing 



oi tne c lusters, models tnat include pure stellar contmuum, are madequate tor age-dat ir 
studies ( Papaderos et al. 1998 : Anders &: Fritze-v. Alvensleben 2003 : Reines et al.lboiol ). 



Some of the causes analyzed in the literature that have been proposed to explain this 
excess are : i) there may be an irnportant contribution from nebular con tinuum and line emis- 



sion ( e.g.. 



le.g- 



Vanzi et al. 



2000 



Papaderos et al.lll998l: IVanzi et al.ll2002l : iReines et al.ll2010l ). ii) hot dust emission 



Reines et al.l l2008l ). iii) the presence of young Red Super Giant 



(RSG) stars in older clust ers, not properly modeled at low metalicity (IMaeder fc Meynet 



2001 



Vazquez et al.l l2007l . a nd references therei n), and iv) other sources as for example 
extended red emission (ERE; IWitt &: Vijh II2004I ). produced by p hotoluminescence pr ocess 
and the presence of a population of young stellar objects (YSOs; lAdamo et al.l l201ll ) still 
surrounded by circumstellar disks. 

Since there is no degeneracy in the color space between age and extinction at young 



ages (Figure 



P, at least qualitatively, model II explains the colors of young star clusters 



(iReines et al. 



20101 ) given that only this model is able to reproduce the trend seen in the 



data, in Figure HI at young ages. 



The presence of nebular continuum and emission lines in the near-IR can have a large 
impact on the inferred properties (affecting the mass determination) of the star clusters in 
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Hll/BCD galaxies, and models that include this effect are the most appropriate in the study 
of young stellar populatioi i with ages 5,6 My r. In fact, near-IR spectra of Tol 35, Tol 3, 
and UM 462 presented by IVanzi et al.l (|2002[ ) clearly show recombination lines of Hel and 
HI (Br7, Bt6 and Pa/3), and excited lines. The line emission detected by Vanzi et al. is 
sufficient to produce their observed broad-band excess. Although hot dust and the presence 
of RSGs may be important factors that can contribute to the near-IR excess, we assume that 
the excess in the SF regions in our sample of Hll galaxies is mainly produced by nebular 
continuum and emission lines, hence we can use our near-IR photometric bands and model 
II, that include this contribution, in order to estimate the properties of the detected star 
clusters/complexes. 



3.1.3. Uncertainties in the age determination 

In order to estimate the uncertainties in our age calculation procedure, we calculate the 
best fit ages while varying the colors by their typical observational errors (as listed in Table 
[2]). In Figure |6] we show the ages obtained from model II for three different extinctions E(B- 
V)=0.0, 0.5 and 1.0 as a function of the mean age (<age>) obtained varying the colors by 
their errors. We showed that in the worst case (if all the errors conspire), they can modify 
the solutions, but not changing significantly the properties of the sample. In this figure, 
we can see that for high extinction the clusters/complexes are better fitted than for lower 
extinctions, producing that the best fits are obtained at higher values of reddening. This is 
an indication that the starburst regions or complexes are dominated by highly extinguished 
and very young star clusters. 

A major concern in our method stems from the fact that we are using near-IR bands 
to obtain the properties of the star cluster/complexes. The determination of absolute 
ages/masses of star clusters and star cluster complexes using near-IR colors, alone, could 
be highly uncertain and highly model dependent. It seems that better constraints of the 
properties of the star clusters/complexes are obtained with a large broad band coverage plus 
emission line equivalent widths. However, we depend on the improvement of population syn- 
thesis models and the contraints of the input parameters and ingredients, as stellar tracks, 
libraries at different metallicities, and inclusion of the effects of nebular emission, dust, etc. 
Given these uncertainties, we cannot obtain absolute ages, but our analysis, based on these 
observations give us a quahtative view of a rather homogeneous star cluster population. Rel- 
ative ages and masses can, however, be useful in order to indicate the mode of SF at galactic 
scales. 

We note that future high-resolution ground-based or space telescope observations using 
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UV-UBVRIJHK broad bands are needed in order to obtain the SED of each star cluster 
or/and complex in our sample of galaxies and then constrain absolute ages and quantify the 
near-IR excess produced by other mechanism, e.g., hot dust, ERE, YSOs, RSGs in metal 
poor ambients, etc. 



3.1. 4- Summary of the obtained properties 

In summary, we found that the star cluster population in Mrk 36 is massive with es- 
timated masses of ~1O'^~^M0. We detected, in this galaxy, a few clusters with masses of 
^IO^Mq. Given our detection limits, lower mass clusters are likely not to be detectable. 
Meanwhile, the star clusters complexes in UM 408 and UM 461 have masses from ~10^Mq 
to ~1O^M0. The age distribution shows that the detected star clusters/complexes are very 
young with ages less than ~10 Myr in Mrk 36 and ~5 Myr in UM 408 and UM 461, respec- 
tively. The two star clusters with colors consistent with ages >10-100 Myr, in Mrk 36, are 
likely old star clusters as is, likely, the case of the cluster #27. Other possibility, is that these 
objects are field galaxies not properly resolved. As we mentioned previously, some of the 
star cluster complexes are highly extinguished and young. In Mrk 36, the less extinguished 
clusters, in most cases, are the youngest star clusters (mainly in complex I). This feature is 
less clear in UM 408 and UM 461, where our resulting properties are averaged over a wider 



region . In the case of UM 408, our results agree with the ones obtained by iLagos et al. 



( l2009l ). where using GMOS-IFU spectroscopy they found that the highest values of extinc- 
tion in the c(H/3) map are displaced from the peaks of Ha emission. The position of the 
peaks of extinction found by Lagos et al. are correlated with the position of our detected 
star cluster complexes. This feature suggests that the current starburst episode is sweeping 
the gas and dust out of the center into the surrounding regions. 

Finally, from a comparison of our results with the ages and masses of other star clusters 
in Hll/BCD galaxies studied in the literature, we see that the ranges are similar. For 
example, the analysis of the cluster population in UM 462 has revealed ages between 4.7 



Adamo et al. 



and 10 Myr (Vanzil 20031 ) with masses range from 1.2 to 7.2xlO^M0. In the case of Haro 11 



(120 lOh found that 30% of the clusters have masses >1O^M0, arguing that these 



clusters qualify as SSCs. Using HST Ima ging Spectrograph (STIS) long-slit far- and near- 



ultraviolet spectra, IChandar et al.l (120041 ) studied a local sample of SSCs in WR starburst 
galaxies, including Mrk 36. They estimated an age <1 Myr for the SSC that coincides with 
the position of our cluster #1, the most luminous cluster in Mrk 36. Given that WR features 
has been reported for Mrk 36 and UM 461 in the literature and since long slit spectra are 
taken centered on the brightest regions, we expect the brightest clusters in Mrk 36 to be 
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really young clusters with upper age limit in the range <2-5 Myr (iGarcia- Vargas et al.lll995l ) . 
If the detected clusters, in Mrk 36, are really young, and given that their ages, masses and 
sizes (with diameters less than 25 pc) are similar to the properties obtained in other star 
clusters, particularly in interacting and star forming galaxies, we suppose that some of the 
most massive star clusters, found in Mrk 36, can be considered candidates to be SSCs. The 
sizes of the regions in the other two H ll galaxies studied here are characteristic of star cluster 
complexes. High spatial resolution observations in space are needed in order to resolve the 
elementary entities that constitute the starburst population in UM 408 and UM 461. 



3.1.5. The spatial distribution of the star cluster / complexes 



In a self-propagating SF model (iGerola fc Seiden 1119781 ) the gas expansion caused by 



stellar-wind and SNe shock waves will triger the next generation of stars. If this scenario of 
SF is plausible, we must derive an age trend of the star clusters, or age sequence of some 
groups of star clusters with respect to the spatial position. To illustrate this, we show in 
Figure [7] the position of the clusters in the galaxy Mrk 36. The star clusters, in this figure, 
are divided into two groups with masses ~1O^M0 and ^IC^Mo, considering three ranges in 
age: 1-5, 5-10 and >10-100 Myr. The spatial distribution of the star clusters/complexes, 
their relative ages and spatial position in the host galaxy show that there is no clear age 
trend, meaning that the young regions were not triggered by the action of the older ones. 

In order to estimate the propagating timescale of the feedback from the star-formation 
activities, we calculate the crossing time tcross=R-/vprop, where R is the size of the system 
and Yprop the propagation or expansion velo city. Typically, this velo city varies from ~10 



kms ^ to ~100 kms ^ in dwarf galaxies (e.g.. Ivan Eymeren et al.l 120071 ). Denoting by R the 



distance between the central or brightest cluster or complex to the other clusters, we found 
that the propagation of SF is possible within the complexes in Mrk 36, for high expansion 
velocities, given that the tcross < Aage (with a mean tcross and Aage~3 Myr in Mrk 36), 
but not on galactic scales because icross > Aage. For low velocities the difference between 
the ages of the central region (in complex I) with respect to the more distant clusters (in 
complex II) are Aage<Ctcross (with a mean tcross ~30 Myr in Mrk 36), suggesting that the 
propagation of SF between the complexes is improbable. 

Given that the star clusters/complexes are practically coeval, the star formation was 
produced simultaneously within time scales of the order of Aage on galactic scales. More 
likely a global SF mechanism is responsible for the present SF activity in galactic scales in 
some Hll galaxies, whereas, self propagating SF on scales of <100 pc is still possible within 
the individual complexes. 



- 16 - 



3.2. Properties of the LSB component 

3.2.1. Surface Brightness Profiles and Structural parameters 



In order to obtain tlie surface brigiitness profile, for eacli galaxy, we fitted ellipses to 
the isophotes using the IRAF task ellipse. First, we manually masked the bright regions 
identified in the Kp-band images, except the central ones, and replacing the values by the 
average of the adjacent regions. Then, we approximated the initial ellipse centers, ellipticities, 
and position angles, allowing these parameters to vary freely with radius during the fitting 
process. When the routine cannot proceed in the iterations as it reaches the lowest surface 
brightnesses, the ellipse task stops, and we fix the center, ellipticity, and the position angle 
with these values to produce our light profiles. 

Figure [S] shows the surface brightness profiles (upper panels) of the isophotal mean 
intensity of each galaxy in the sample, considering pixels with la above the background, for 
the J, H, Kp and Br7 filters. Br7 is in arbitrary units. We measure the fainter level limits to 
be /ij ~23.50, 23.10 and 23.50 mag arcsec"^, ^23.00, 23.10 and 23.40 mag arcsec"^, and 
^Kj, ^22.30, 22.50 and 23.00 mag arcsec'^ for Mrk 36, UM 408, and UM 461, respectively. 
The J-H and H-Kp color profiles are shown in the lower panels of Figure [Hi We note that the 
color profiles show relatively constant underlying values at intermediate radii with a color 
gradient at large radii, indicating that their stellar populations must be fairly homogeneous 
in the main body of the galaxy with the presence of a old stellar population at large radii. 
In Figure IH] we show the pixel-to-pixel color map of the galaxies in order to compare the 
color spatial distribution of these galaxies with the profiles previously obtained. Again, we 
note that the color spatial distributions are relatively constant through the main body of the 
galaxies. The surface photometry of the galaxi es were corrected on ly for Galactic extinction 
using the relation with Ri/ =Av/E(B-V)=3.1 (jCardelli et al.lll989l ) and adopting the values 
from the extinction map of Schlegel et al. ( IQQsl ) showed in Table [T] Column (7). We did not 
apply any smoothing procedure to our data in order to obtain the surface brightness at large 
radii. 



We see that the out er part of the light profiles in Figure [8] are wel l represented by an 
exponential model (e.g., iPapaderos et al.l Il996l : iTelles fc TerlevichI llQQTt ICairos et al.l l2003l : 
Noeske et al.ll2003l . |2005| ) thus other functions are not necessary. We can express this profile 
in terms of the surface brightness, so we can find that yu(r)=/io^A+(l-086/a;A)r, where /io,A 
corresponds to the central surface brightness and ao,{j,j/,Kp) the scale length with A=J, H 
and Kp. Column (1) of Table S] gives the name of the galaxies. Columns (2), (3) and (4) give 
the /io,(j,H,Kp)! C(o,iJ,H,Kp) parameters and raLSB,{j,H,Kp) magnitudes of the LSB component 
described below in § I3.2.2[ respectively. These parameters will be used in § 13.2.21 in order 
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to calculate the integrated properties of the LSB component. 



3.2.2. The integrated properties of the LSB component 



The total apparent magnitude of the LSB component can be obtained using that m.LSB 
= fiQ - 5 log(Q;) - 1.995 - 2.5 log(l-e), where e correspond to the ellipticity derived from the 
surface brightness profile (with constant values of £=0.40, 0.30 and 0.33 for the galaxies Mrk 
36, UM 408 and UM 461, respectively). So, we calculated the J-H and H-K colors of the 
LSB component from the results of the exponential fits. The colors of the LSB compori ent in 
Mrk 36 are similar to those obtained previously in the literature bvlCairos et al. (|2003[) with 
J-H=0.37 and H-K,=0.28 from exponential fits. For UM 461, iNoeske et all (]2003|) found 
J-H and H-K^ colors, of the LSB component, that disagree and agree with our values at 
la level, respectively. The reason for this disagreement could be the uncertainties in the 
sky estimation and/ or the aperture diffe rences in both studies causing the surface brightness 
profile obtained by iNoeske et al.l ( l2003l ) to be deeper with the radius compared with our 
profiles, thus producing higher slopes of the profiles at large radii. In summary, the LSB 
component of our sample of galaxies show blue J-H colors and red H-K that suggest a 
photometrically dominant stellar populatio i i of ag es > 10^ yr. Our interpreta tion agrees with 
the age estimates given by iRaimann et al.l (120001 ) and IWestera et al.l (120041 ) for a sample of 
Hll galaxies (in which they include UM 408 and UM 461 in their analysis) by means of 
spectral population synthesis. 

Finally, we estimated the stellar mass (M^,) of the LSB coni ponent for each galaxy using 
the M/L relationship for the H-band (see iLopez-SanchezI |2010| ). We found that the stellar 
mass for our sample of galaxies are log(M^)= 7.89, 8.33, and 7.88 Mq for Mrk 36, UM 408, 
and UM 461, respectively. These stellar masses are typical for BCD galaxies calculated using 
different photometric bands, and dominated by the contribution of the intermediate to old 
stellar population (e.g., 3xlO^M0 for Mrk 36: lAmorm et al.ll2009l ). 



Discussion 



see 



Star formation o ccurs when the local molecular gas density exceeds a certain threshold 



Leroy et al.ll2008l . and references therein), given first possibly by collisions of gas clouds. 



due to turbulence or gravity, and then triggered by the action of massive star evolution, e.g. 
stellar winds from star clusters and SNe. In many cases, the current star-formation may be 
triggered by external agents, as in the case of the tidal forces by a neighboring galaxy or 
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mergers. However, for a fraction of isolated, less luminous and compact Hll galaxies, this 
star-formation activity may occur solely by internal processes. 

The SF activity r nay form a fract ion of bound star clusters that will survive their infant 
mortality and evolve. iBastian I (120081 ) has defined the value of the present cluster formation 
efficiency F in a host galaxy as CFR/SFR, where CFR = Mtot/At is the present cluster 
formation rate at which the galaxy produces a total cluster mass Mtot in a given age interval 
At. On the other hand, iGoddard et al.l (120101 ) found a correlation between the value of the 
present cluster formation efficiency and the SF density of the host galaxy in a sample of 
starburst galaxies, T{%)=29J2'gp%i^ Mq jt~^ kpc~^, where T^sfr is the total SFR per unit of 
area. 

The first step in the determination of the star cluster formation efficiency in our sample 
of g alaxies is the dete rmination of the total SF rate (SFR). We calculate the total SFr[^ using 
the iKennicuttI (119981 ) relationship for L(Ha). We use the integrated B.(3 flux ob tained using 



contin uum-free emission line broad band images, with L(Ha)/L(H/3)=2.87, by iLagos et al. 



(120071) for Mrk 36 and UM 461, and integrated integral fleld unit Ha flux obtained by 
Lagos et all J2009h for UM 408. So we flnd that the total SFR(Ha)~ 0.078 Mq yr'^ in Mrk 
36, 0.017 Mq yr"^ in UM 408, and 0.085 Mq yr'^ in UM 461. The SFR in our sample 
of galaxies is very low compared with other dwarf g alaxies, such as Haro 11 (22 M0 yr~^; 



Adamo et al.ll2010h and NGC 1519 (0.3626 M© yr'^; IGoddard et al.ll2010f ) and in fact with 



respect to more irregular and lumino us Hll/BCD galaxies, such as Tol 9 with a SFR(Ha) = 



1.82 Mq yr ^ (ILopez-Sanchea l2010l ). Assuming a starburst area of ~2 kpc^ in Mrk 36, 



~8.16 kpc^ in UM 408 and ~1.62 kpc^ in UM 461, we found that the integrated SFR per 
unit of area in Mrk 36 is Esfr=0-039 Mq yr^^ kpc^^, in UM 408 is Esfr=0-00'2 Mq yr'^ 
kpc~^ and in UM 461 is £5^7/2=0. 052 Mq yr^^ kpc^^, r espectively. We can n ow calculate 
the expected value of F using the correlation found by IGoddard et al.l (120 lOl ). Thus, we 
obtain that the cluster formation efficiency is ~13% in Mrk 36, ~7% in UM 408 and ~14% 
in UM 461, respectively. Hence, the star cluster efficiency in our low luminosity Hll galaxies 
is approximately ~10%. This curr ent SF efficiency is lower than the ones found in more 
luminous BCD galaxies as Haro 11 ( lAdamo et al.ll2010l ) with an efficiency of ~38%. 



Using the properties of our detected star clusters/complexes showed in Table El we 
calculate a total cluster /complex mass, by summing the masses, of 54.50 xIO^Mq in Mrk 
36 (with a total mass in clusters >10^Mq of 44.69x10^ Mq), 13.81 xIO^Mq in UM 408 and 



^We multiplied the SFR by a factor 0.67 obtained from the comparison between the SFR of our galaxies 
using the IKennicuttI (|l998l ) and Calzetti et all (|2007l ) relationships This due to differences in the stellar IMF 
assumptions given that the IKennicuttI (|l998l) relationship is based on a Salpeter IMF and the models which 
we adopted are based on a Kroupa IMF 
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43.09 xlO^M© in UM 461. So, we calculate that the fraction of the total cluster mass with 
respect to the LSB host galaxy mass in Mrk 36, considering our complete range in ages, is 
equal to -0.007, ~0.006 in UM 408 and ~0.006 in UM 461. This implies that the inferred 
total current SF ma ss is of the ord e r of 5 ,1% of the underlying galaxy mass, in agreement 
with the estimate of IWestera et al.l (120041 ) that the past history of SF in the galaxies were 
more active than the present one. Additionally, we calculate the SFR per unit of area for the 
individual star clusters/complexes using our measured Br7 emission. From these, we obtain 
that in the brightest regions (clusters t^I and #5 in Mrk 36 and complex #5 in UM 461) the 
SFR per unit of area i s compatible with the values observed in starburst galaxies with values 
>0.1 M0 yr"^ kpc~^ (IBastian et al. II2005I . and references therein). While, the majority of 
the star clusters/complexes have values of the order of —0.01 Mq yr~^ kpc~^. If we use only 
the total mass in clusters >10'^Mq to calculate SF efficiency, we obtain that in Mrk 36 for 
an age of 20 Myr our result agre e , with in the uncertainties, with the one obtained using 
the relationship of iGoddard et al.l (120101 ). However, in the case of UM 408 and UM 461 the 
value of r is extremely high, indicating that in these cases the star forming complexes are 
not resolved into individual clusters, resulting in a overestimation of the F parameter. 

We can also calculate the gas consumption time scale (Tgas) or how long it would take 
before all the gas in the galaxies will be consumed at the current SFR. The consumption time 
scale is defined as the ratio between the available gas and the current SFR , r^ns = M^n, , /SFR . 



gas 



Mni+MHe+MH^ ^ 2x M^r (ILerov et al. 



Assuming the total amount of gas M 
we obtain that r„.. - .5 Gyr in Mrk 36 (with Mi^/=2.OxlO'^M0; iThuan fc Martin 



Bravo- Alfaro et al.l 120041 ) . Tgas ~ 77 Gyr in UM 408 (with Mf/f =6.53xlO^M^. ijai^t;r tji ai.i 
2002h and Tgas ~ 2 Gyr in UM 461 (with Mh/=O.98x1OSM0; ISmoker et aPboOol ). The 



2005h 



1981 



Salzer et al. 



consumption timescales in Mrk 36 and UM 461 are signi ficantly less than a Hub ble time and 
comparable with the times observed in spiral galaxies ( iKennicutt et al. Ill994j ). So the SF 
cannot be sustained for the entire history of the galaxies, which indicates that these objects 
undergo a few or several short bursts of SF. While UM 408 has a gas consumption timescale 
longer than a Hubble time, indicating that likely the SF is relatively constant through the 
history of the galaxy or the HI halo may not be spatially available for the current SF. 

Our findings also seem to indicate that the SF mode in our sample of low luminosity 
Hll galaxies is clumpy, similar to other dwarf galaxies (e.g., NGC 1569 and SBS 0335-052). 
These complexes or star-forming knots a re formed by a few m assive star clusters with masses 
>1O'^M0 and high SFR per unit of area. iMelena et al.l (120091 ) found that the observed trends 
in the number and mass of the SF regions, in a sample of dwarf galaxies, is independent of 
the local environment and even the surround ing galaxy mass and i s given by a mass fu nction 
which stochastically favored SF in clusters (lAdamo et al.ll2010l ). iBillett et al.l ( 120021 ) show 
that SSCs may require special (or fortunate) circumstances to form in dwarf galaxies, but 
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when they do, they are very massive (>1O^~^M0) and form clumps or groups of similar 
ages. This suggests that these clumps are likely formed in localized regions of high pressure 
triggered by large scale ambie nt gravitational instabilitie s, given that in dwarf galaxies most 
of the ISM is at low pressure (lElmegreen fc Hunterll2000l ). The lack of external perturbers in 
the most compact and isolated galaxies indicates that an additional mechanism other than 
tidal interactions must be considered to explain this current SF activity. This mechanism 
may be related to the overall physical conditions of the ISM, particularly the gas surface 
densities, in conjunction with stochastic effects, that allow SF to take place. Alternatively, 
the low current SFR implies that a burst or a triggering is not necessary, simply that the SFR 
has been relatively constant. However, tidal interactions or mergers are likely the primary 
agent to trigger the cu rrent SF in luminous and more disturbed H ll galaxies as suggested by 
their morphology (e.g.. lTelles fc Terlevichlll995l : iLagos et al.ll2007l : iLopez-Sanchez &: Esteban 



20081 ). 



5. Conclusions 

In this paper, a sample of three Hll galaxies (Mrk 36, UM 408, and UM 461) has 
been analyzed in order to study their stellar populations (star cluster complexes and the 
underlying host galaxy or LSB component) using new near-IR high spatial resolution images 
obtained on the Gemini North telescope. In our analysis we used models that include the 
contribution of stellar continuum, nebular continuum and emission lines. Our conclusions 
can be summarized as follows: 

1. The presence of nebular continuum and emission lines in the near-IR produces an 
excess in the observed SED in young star cluster/complexes. This excess, can have a 
large impact in the inferred properties of the star clusters in Hll/BCD galaxies and 
models that include this effect are the most appropriate in the study of young stellar 
population with ages <6 Myr. 

2. We found that the star cluster population in Mrk 36 shows masses of ~1O^~^M0 with a 
few detected star clusters with masses of ~1O^M0 distributed in the main body of the 
galaxy. The star cluster complexes in UM 408 and UM 461 have masses from ^IO^Mq 
to ~1O^M0. The age distribution shows that the detected star clusters/complexes are 
very young with ages of a few Myr. Two likely old star clusters with colors consistent 
with ages >10-100 Myr have been detected in Mrk 36. The fraction of recent SF in 
bound clusters and/or complexes more massive than 1O^M0 is compared to the current 
total SFR is about 10% in our sample of galaxies. 
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3. The spatial distribution and ages of the star cluster/complex population seems to in- 
dicate that SF is clumpy and simultaneous. We propose that the current SF activity 
in our sample of low luminosity H II galaxies is triggered by some internal mechanism 
instead of tidal interactions. This mechanism of SF may be related to the overall phys- 
ical conditions of the ISM that produce the increase of surface densities in conjunction 
with stochastic effects within a time scale comparable to the mean age differences of 
the massive star cluster complexes. 

4. The LSB component of our sample of galaxies have near-IR colors representative of 
evolved stellar population of at least >10^ yr. We found that the stellar mass of this 
component for our sample of galaxies arc Iog(M*)= 7.89, 8.33, and 7.88 Mq for Mrk 36, 
UM 408, and UM 461, respectively. The fraction of the total cluster mass with respect 
to the LSB hosting galaxy mass in our sample of galaxies, considering our complete 
range in ages, is less than 1%. 
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Fig. 2. — Kp and Br7 images (not continuum subtracted) of the galaxies Mrk 36 (top), UM 
408 (middle) and UM 461 (bottom) with NIRI at Gemini North. The field of view of each 
image is 48"x48". The images are displayed on a logarithmic scale. The arrows in the Kp 
image of Mrk 36 point to three background galaxies close to the galaxy. North is at the top 
and east to the left. 
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Fig. 3.— Images of Mrk 36, UM 408, and UM 461 in the Kp filter. Circles mark the 
position (and the apertures used for photometry) of each star cluster/complex identified in 
the galaxies. The clusters have been labeled with the designation used throughout the paper 
(see Table [3]). In the upper panel of Mrk 36 we show the position of each cluster identified 
in the northern region of the galaxy. The circle at th e corner of this panel represents the size 
of the nebular region of 30 Doradus (jWalbornlll99ll ). The lower panel shows the clusters in 
the central region of Mrk 36 (dashed square in the upper panel). 
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Fig. 4. — Observed near-IR color-color diagram J-H vs. H-K for the star clusters/complexes 
in our sample not corrected for extinction. The lines show the evolutionary tracks of these 
colors from STARBURST99 (black line; this model includes stellar and nebular continuum) 
and GALEV models (orange line; this model includes stellar continuum, nebular continuum 
and the contribution of nebular emission lines) for a metallicity Z=0.004. Red triangles 
represent the observed values for clusters in Mrk 36, the blue stars represent the values 
for UM 408 and green circles represent the values for clusters/complexes in UM461. Filled 
symbols indicate the detection of Br7 in the regions. The error bars show the average error 
value for the colors. Open circles along the tracks indicate ages of 1, 3, 6, 10, and 100 Myr. 
An additional circle indicating an age of 1 Gyr is included in the model of STARBURST99. 
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Fig. 5. — Extinction, age and mass distribution of our sample of star clusters/complexes, 
obtained using the STARBURST99 (model I; black distribution) and GALEV (model II; 
orange distribution) models. For more details see § I3.1.2I 
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Fig. 6.— Best fit ages (listed in table [2]) from model II for Mrk 36 (red), UM 408 (blue) 
and UM 461 (green) for different extinctions E(B-V) as a function of the mean age obtained 
varying the colors by their typical observational errors. The errors in this figure were calcu- 
lated from the difference between the ages from the observed colors and the ones obtained 
by displacing these colors by their la uncertainties. 



- 32 - 



1-5 5-10 >10-100 Myr 



10* M„0 O 
10^ M„ O 



15 



complex II 



o 



o 



10 



o 



d o 




o 



o 



5 



complex I 



252 pc 

I I I \ I I I I \ I I I I \ I I I L 

5 10 15 

arcsec 



Fig. 7. — Spatial distribution of the star clusters in Mrk 36. The different circle sizes 
correspond to masses ^lO^M© and ^IO^Mq considering three ranges in age: 1-5, 5-10 and 
>10-100 Myr. 
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Fig. 8. — Upper panels: surface brightness profiles of Mrk 36, UM 408, and UM 461 in J, H, 
and Kp (corrected for galactic extinction) and Br7 in arbitrary units. We considered pixels 
with 1(7 above the background. For a better visualization, the H and Kp profiles are shifted 
by -1 and -2 mag, respectively. Bottom panels: J-H and H-Kp color profiles. 
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Fig. 9. — J-H and H-K color maps of the galaxies Mrk 36 (top), UM 408 (middle), and UM 
461 (bottom). Black contours in the J-H maps corresponds with J-band contour maps and 
in the H-K maps corresponds with Kp-band contour maps. Orange contour are from Br7 
images. 



Table 1. Sample data and NIRI observations. 



Object a 5 Vel. 12+log(0/H) c(H/3) E(B-V)Ga! Date of Filter Exposure Air mass Ca 

(J2000) (J2000) (km/s) Observation Time (s) (mag) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



Mrk 36 11:04:44.0 +29:07:48 646 7.81'' 0.02=' 0.031 2005-11-24 



UM 408 02:11:23.4 -1-02:20:30 



UM 461 11:51:33.1 -02:22:22 



(12.95 Mpc) 








J 


6x(3x60) 


1.37 


23.93±0.05 










H 


6x(4x65) 


1.27 


24.01±0.05 










Kp 


12x(3x75) 


1.15 


23.57±0.06 










Br7 


6x(3x40) 


1.08 




3598 


7.87'' 


0.93'' 0.037 


2005-08-02 










(45.70 Mpc) 








J 


6x(3x60) 


1.26 


23.91±0.05 










H 


6x(4x65) 


1.19 


23.97±0.03 










Kp 


12x(3x75) 


1.10 


23.49±0.04 










Br7 


6x(3x40) 


1.06 




1039 


7.78'' 


0.12" 0.018 


2005-12-29 










(19.20 Mpc) 








J 


6x(3x60) 


1.62 


23.95±0.07 










H 


6x(4x65) 


1.46 


24.05±0.05 










Kp 


12x(3x75) 


1.28 


23.54±0.02 










Br7 


6x(3x40) 


1.18 





00 



' llzotov fc ThuaiJ l|l998h 
' jLagos et al ] l|2009l l 



Note. — Column (1) galaxy name. Columns (2) and (3) a and S coordinates (J2000), respectively. Column (4) heliocentric velocity (vel.) and the 3K 
CMB corrected distance from NED. Columns (5) and (6) ox ygen abundanc e and extinction adopted in this work. Column (7) date of observation. Column 
(8) Galactic extinction obtained from the extinction maps of lSchlegel et al.l lll998l V Column (9) filter used in each observation and finally columns (10) and 
(11) exposure time in s and the mean air mass of each observation, respectively. Column (12) shows the instrumental zero points C^ with A= J, H and Kp 
obtained for each run of observation. 
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Table 2: Integrated magnitudes of our sample of galaxies. 



Object 


J 


H 


Kp 




(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


Mrk 36 


14.46±0.05 


14.23±0.05 


13.70±0.06 


UM 408 


15.94±0.07 


15.72±0.06 


15.28±0.05 


UM 461 


15.04±0.09 


14.73±0.07 


14.53±0.05 



3. Observed aperture photometry in the near-IR bands J, H and Kp and measured properties of the stellar 

clusters. 



Name Cluster Aperture radii J H Kp E(B-V) Age Mass Log L(Br7) 
(arcsec) (mag) (mag) (mag) (mag) (Myr) (xlO^M©) (erg s~^) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 

Mrk 36 



It 


0.44 


18.64±0.05 


18.58±0.05 


18.09±0.06 


0.00 


-0.40 


1.99 - 


2.8 


5.46 - 


2.53 


36 


.256 


2 


0.39 


20.01±0.05 


20.13±0.05 


19.44±0.06 


0.00 


- 0.35 


1.25 - 


1.1 


1.73 - 


0.64 






3t 


0.28 


19.98±0.05 


19.58±0.05 


19.40±0.06 


0.50 


- 0.75 


10.11- 


8.6 


4.08 - 


3.62 


35, 


.268 


4t 


0.30 


19.88±0.05 


19.51±0.05 


19.24±0.06 


1.00 


- 1.00 


5.07 - 


5.5 


3.49 - 


1.77 


35 


.818 


5t 


0.29 


19.99±0.05 


19.68±0.05 


19.46it0.06 


0.95 


- 0.75 


6.09 - 


5.5 


3.56 - 


1.33 


35 


.912 


6t 


0.37 


19.84±0.05 


19.45±0.05 


19.19±0.06 


1.00 


- 1.00 


5.07 - 


5.5 


3.66 - 


1.86 


35 


.314 


7t 


0.29 


19.87±0.05 


19.52±0.05 


19.23±0.06 


1.00 


- 0.90 


4.89 - 


5.0 


3.81 - 


1.64 


35 


.538 


8t 


0.40 


18.83±0.05 


18.62±0.05 


18.45±0.06 


0.60 


- 0.50 


5.62 - 


5.9 


6.78 - 


3.60 


35 


.072 


9 


0.24 


20.74±0.05 


20.54±0.05 


20.21±0.06 


0.20 


- 0.65 


2.95 - 


4.2 


0.97 - 


0.62 






10 


0.28 


19.99±0.05 


19.61±0.05 


19.57it0.06 


0.00 


- 0.00 


~1000 - 


- 19.3 


45.62 - 


- 4.50 


35, 


.523 


11 


0.29 


19.74±0.05 


19.62±0.05 


19.50±0.06 


0.30 


-0.25 


5.43 - 


6.5 


2.18 - 


1.70 






12t 


0.24 


20.48±0.05 


20.28±0.05 


20.19±0.06 


0.20 


- 0.25 


7.85 - 


8.3 


1.18 - 


1.44 


34, 


.273 


13 


0.29 


20.57±0.05 


20.43±0.05 


20.06±0.06 


0.00 


-0.45 


2.78 - 


4.0 


0.99 - 


0.70 






14 


0.32 


20.62±0.05 


20.37±0.05 


20.05±0.06 


0.45 


- 0.75 


3.67 - 


4.3 


1.37 - 


0.71 






15 


0.43 


20.97±0.05 


20.83±0.05 


20.21±0.06 


0.30 


- 0.75 


1.25 - 


1.0 


0.91 - 


0.36 






16 


0..38 


21.10±0.05 


20.97±0.05 


20.50±0.06 


0.00 


- 0.65 


1.64 - 


3.5 


0.60 - 


0.42 






17t 


0.56 


19.39±0.05 


19.10±0.05 


18.82±0.06 


0.75 


- 0.80 


4.46 - 


4.8 


5.52 - 


2.26 


35, 


.374 


18 


0.52 


19.59±0.05 


19.54±0.05 


18.87±0.06 


0.20 


- 0.65 


1.25 - 


1.0 


3.08 - 


1.21 






19t 


0.44 


19.80±0.05 


19.39±0.05 


19.21±0.06 


0.50 


- 0.80 


10.23 - 


- 8.5 


4.77- 


4.40 


33, 


.516 


20t 


0.40 


19.99±0.05 


19.70±0.05 


19.39±0.06 


0.70 


- 0.70 


4.16 - 


4.5 


3.03 - 


1.33 


35 


.497 


21 


0.40 


20.40±0.05 


19.96±0.05 


19.80±0.06 


0.25 


- 0.55 


26.92- 


9.6 


3.39 - 


2.26 






22 


0.27 


21.19±0.05 


21.05±0.05 


20.76±0.06 


0.00 


- 0.50 


2.95 - 


4.3 


0.54 - 


0.34 






23 


0.41 


20.28±0.05 


19.92±0.05 


19.61±0.06 


1.00 


- 1.00 


4.57 - 


4.8 


2.90 - 


1.18 






24t 


0.47 


19.75±0.05 


19.35±0.05 


19.14±0.06 


0.85 


- 0.95 


15.85 - 


- 7.9 


9.72 - 


4.79 


35, 


.060 


25t 


0.36 


20.68±0.05 


20.40±0.05 


20.13±0.06 


0.75 


- 0.75 


4.57- 


4.8 


1.64 - 


0.66 


34, 


.859 


26t 


0.33 


20.75±0.05 


20.50±0.05 


20.10±0.06 


0.05 


- 0.85 


1.25 - 


4.1 


0.89 - 


0.77 


34, 


.829 


27 


0.33 


20.89±0.05 


20.54±0.05 


20.61±0.06 


0.00 


- 0.00 


~1000 - 


- 17.9 


17.20- 


1.61 






28t 


0.36 


20.39±0.05 


20.07±0.05 


19.74±0.06 


0.80 


- 0.95 


4.21 - 


4.5 


2.31 - 


0.98 


35, 


.261 


29t 


0.40 


19.99±0.05 


19.76±0.05 


19.32±0.06 


0.10 


- 0.80 


1.25 - 


4.0 


1.87 - 


1.60 


33 


.584 


30 


0.41 


20.52±0.05 


20.30±0.05 


19.93±0.06 


0.00 


-0.70 


2.39 - 


4.1 


0.99 - 


0.85 






3lt 


0.25 


20.11±0.05 


19.92±0.05 


19.66±0.06 


0.25 


- 0.55 


3.84 - 


4.6 


1.74 - 


0.92 


34 


.916 


32 


0.50 


20.57±0.05 


20.31±0.05 


19.93±0.06 


0.35 


- 0.85 


2.88 - 


4.2 


1.32 - 


0.86 







Table 3 — Continued 



N&mG 




ADprtnrp Tarlii 


J 


H 




E(B-V) 


Ap'p 


Mass 














( vn Q rr I 


(mag) 


^'^Avr^ 

(LVLyr) 




(.erg s ; 










C5") 


("J 


(1 ) 






no") 




33+ 


0.26 


20.14±0.05 


20.06±0.05 


19.71±0.06 


0.00 - 0.35 


2.88 - 4.0 


1.44 - 0.95 


34.997 


UM 408 






















1 


0.40 


20.52dz0.05 


20.05±0.03 


19.74±0.04 


1.00 - 1.00 


9.66 - 5.1 


43.69 - 13.48 






2t 


0.53 


18.95±0.05 


18.53±0.03 


18.24±0.04 


1.00 - 1.00 


5.07 - 5.2 


110.22 - 54.26 


36.087 




3t 


0.27 


20.49±0.05 


20.11±0.03 


19.82±0.04 


1.00 - 1.00 


4.78 - 5.1 


28.63 - 12.56 


35.457 




4t 


0.33 


20.19±0.05 


19.88±0.03 


19.44±0.04 


0.25 - 1.00 


1.25 - 4.1 


21.88 - 18.57 


35.495 




5t 


0.40 


19.77±0.05 


19.40±0.03 


18.99±0.04 


0.90 - 1.00 


3.67 - 4.2 


53.64 - 27.05 


35.762 




6 


0.17 


21.63±0.05 


21.18±0.03 


20.92±0.04 


1.00 - 1.00 


7.94 - 8.2 


9.81 - 12.21 




UM 461 






















1 


0.66 


19.87±0.07 


19.63±0.05 


19.37±0.02 


0.60 - 0.70 


4.4 - 4.7 


6.94- 2.85 






2t 


0.86 


17.34±0.07 


17.50±0.05 


16.86±0.02 


0.00 - 0.15 


2.0 - 1.0 


38.74 - 14.33 


37.206 




3 


0.42 


20.00±0.07 


19.85±0.05 


19.39±0.02 


0.00 - 0.65 


1.3 - 3.7 


3.74 - 2.88 






4 


0.62 


19.88±0.07 


19.44±0.05 


18.85±0.02 


0.85 - 1.00 


1.3 - 2.9 


8.40 - 3.50 






5 


0.55 


19.86±0.07 


19.83±0.05 


19.38±0.02 


0.00 - 0.30 


2.7 - 2.9 


3.96 - 1.65 






6 


0.49 


20.57±0.07 


20.25±0.05 


20.00±0.02 


1.00 - 0.85 


5.5 - 5.4 


3.98 - 1.82 






7 


0.50 


20.27±0.07 


20.08±0.05 


19.82±0.02 


0.30 - 0.60 


3.9 - 4.6 


3.53 - 1.79 






8 


0.57 


20.33±0.07 


20.05±0.05 


19.72±0.02 


0.65 - 0.85 


3.1 - 4.4 


3.33 - 2.16 






9 


0.62 


20.15±0.07 


20.01±0.05 


19.79±0.02 


0.15 - 0.40 


4.1 - 4.7 


3.71 - 1.73 






10 


0.66 


20.14±0.07 


19.81±0.05 


19.52±0.02 


0.90 - 0.90 


4.6 - 4.9 


6.65 - 2.75 






11 


0.61 


20.34±0.07 


20.03±0.05 


19.83±0.02 


0.85 - 0.80 


7.2 - 7.0 


5.32 - 4.10 






12 


0.53 


20.48±0.07 


20.30±0.05 


20.03±0.02 


0.25 - 0.60 


3.8 - 4.5 


2.73 - 1.46 






13 


0.60 


20.60±0.07 


20.20±0.05 


19.93±0.02 


1.00 - 1.00 


5.1 - 5.5 


4.06 - 2.07 





Note. — Column (1) galaxy name. Column (2) identification number of the star clusters/complexes. The ^ symbol indicates that Br7 
emission were measured in the cluster. Column (3) aperture considered to obtain the photometry. Columns (4), (5) and (6) observed 
photometry in magnitude of each star cluster /complex in the filters J, H and Kp, respectively. Columns (7), (8) and (9) extinction E(B-V), 
age in units of Myr and the mass in units of Mq for models I and II, respectively. Finally Colum (10) give the Br7 luminosity. 
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Table 4: Structural parameters from the exponential fits to the host galaxies of Mrk 36, UM 
408 and UM 461. Column (1) galaxy name. Columns (2), (3) and (4) //o,(j,i7,Kj,), Oio^{j,H,Kp) 
parameters and raLSB,{j,H,Kp) magnitudes. 



Object 










(mag arcsec"^) 


(mag arcsec"^) 


(mag arcsec"^) 




ao,j 








(arcsec) 


(arcsec) 


(arcsec) 




T^LSB,J 




T^LSB,Kp 




(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


Mrk 36 


18.64 


18.25 


18.38 




3.79 


3.75 


4.43 




14.31 


13.94 


13.70 


UM 408 


18.90 


18.03 


18.19 




1.91 


1.48 


1.84 




15.89 


15.58 


15.27 


UM 461 


19.51 


19.02 


19.11 




3.69 


3.37 


3.84 




15.11 


14.82 


14.62 



